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Abstract. A Skolem sequence is a linear arrangement of the multiset, {1, 1, 2, 2, . . . , n, n} such that if 
r S [n] appears in positions i and j, then \i — j\ = r. We first translate the problem to a particular set of 
^>f^ ■ perfect matchings, then apply the method of generating trees for open arc diagrams to generate exhaustively 

all Skolem sequences of a given size. Tracking the arc length between pairs of vertices in an arc annotated 
' diagram is the central task. Although we do not surpass previously known enumerative results, this method 

, drastically reduces the search space compared to previously known methods. 

00 '• 1- Introduction 

CN ■ 

Arc diagrams, also known as arc annotated sequences, are structures that encode a variety of combinatorial 
classes, including matchings 0, colored matchings [5], set partitions [H], permutations [S], labelled graphs 
' [S]; tangled diagrams Til and RNA substructures |T7]. Each arc diagram is a row of increasingly labelled 

I vertices from 1 to n with some arcs between them, restricted according to the combinatorial class being 

represented. Much recent work has focused on illustrating the equidistribution of two statistics that arise in 
arc diagrams: fc-crossings (a set of k mutually crossing arcs) and fc-nestings (a set of k mutually nesting arcs) 
( [H] ) [H] ; [S] 1 [H] ; [H] ) ■ Also of interest has been the enumeration of combinatorial classes parameterized 
by these statistics ([5], [B], [SO])- More recently, in [5], a generalized version of an arc diagram was used 
for the first time to construct generating trees and functional equations for fc-nonnesting set partitions and 
permutations. This generalization was the open arc diagram, which allowed for arcs to have left endpoints 
^ . but no right endpoints, and was introduced for bijective purposes in |19| . 

\l ' Here we represent a different combinatorial class, Skolem sequences, with arc diagrams, and utilize open 

. arc diagrams to build a generating tree that exhaustively constructs all Skolem sequences of order n using a 

\^ I strategy similar to that seen in . While this direct generation scheme is restricted in practical effectiveness 

due to memory limitations, we feel this method is valuable for three reasons: 

, (1) This is a new approach that recursively constructs Skolem sequences; 

' (2) The strategy is a systematic approach for various generalizations of Skolem sequences; 

(3) This method fits within the larger picture of open arc diagrams being used as a tool for generation 
and enumeration. 



2. Plan of the paper 

The paper will proceed as follows: in Section [3] Skolem sequences will be defined and background from 
the literature will be put into context. In Section |4] we will depict Skolem sequences as arc diagrams, 
and introduce a generalization of them called open Skolem sequences which will be critical to our generation 
scheme. In Section[S]we present a label for open Skolem sequences that allows us to formalize a succession rule 
for generating all descendants of an open Skolem sequences. We finish with a discussion in Section [S] which 
places memory limitations into context, highlights our contribution to the study of Skolem sequences and 
potential avenues for other generalizations of Skolem like sequences. Lastly Skolem sequences are viewed as a 
subclass of perfect matchings and an observation is made regarding their corresponding oscillating tableaux 
as seen in 0. 

3. Skolem sequences 

Skolem sequences were introduced in 1957 [22], when Thoralf Skolem asked about the possibility of 
partitioning the set {1, 2 . . . , 2n} into n pairs (a^, br) where br — ar ^ r for r = 1, 2, . . . , n. Such a partition 
forms what is now referred to as a Skolem sequence. 
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Definition 1. A Skolem sequence of order n is an integer sequence w — (si, S2, ■ ■ ■ , S2n) of the multiset 
{1, 1, 2, 2, . . . , n, n} of size 2n satisfying the following conditions: 

(1) For all k £ {1,2, .. .n}, there exists exactly two elements Si, sj in w such that Si — sj = k. 

(2) // Si = Sj — k, and i < j , then j ~ i = k. 

We let S denote the set of all Skolem sequences and iS„ to be of order n. 

Example 1. When n = A, w — (3, 4, 2, 3, 2, 4, 1, 1) £ Si is a Skolem sequence. 

Skolem sequences have immediate further applications in design theory, particularly, in Steiner Triple 
Systems. Let STS{n) denote the set of pairs (V,B) where \V\ = v and B consists of 3-subsets ('blocks') of V 
such that any 2-subset of V is included in exactly one block. These are the Steiner Triple Systems of size n. 
In 1958, Skolem [23^ showed how to build an element of STS{Qn-\- 1) given an element of 5„. Specifically, if 
w — {si, . . . S2n) £ Sn, an ST S (fin +1) can be constructed with base blocks of the form {x,x + k,x+j + n), 
where Si — Sj ~ k, i < j and x £ {0, 1,2..., 6n}. 

Example 2. Let n = 4. The Skolem sequence w = (3, 4, 2, 3, 2, 4, 1, 1) £ 1S4 can be written as: 

i= 12345678 
= 34232411 
We build base blocks of a STS(25) from this example using x = 0. Notice that Si — Sii — 3, so our first base 
block is (0,0 + 3,0 + 4 + 4) = (0,3,8). Similarly, S2 = sg = 4 gives the base block (0,0 + 4,0 + 6 + 4) = (0,4, 10); 
S3 = S5 = 2 gives (0, + 2, + 5 + 4) = (0,2,9) and 57 = sg = 1 gives (0, + 1, + 8 + 4) = (0,1,12). Thus 
our base blocks for creating the STS{25) are 

(0, 3, 8), (0, 4, 10), (0, 2, 9), and (0, 1, 12). 
Adding integers 1, 2, . . . , 24 to each element of each base block and reducing mod 25 constructs an STS{25). 

Beyond STSs, Skolem sequences can also be applied to the areas of starters [21], balanced tertiary designs 
[4] , factorizations of complete graphs [T2| and labellings of graphs [2 . For this reason, generation of Skolem 
sequences of arbitrary order is valuable, although most existing algorithms only construct a single, or at least 
limited number, of Skolem sequences of order n. Despite this, since its introduction, most work has centered 
around the existence of Skolem and Skolem-type sequences, the generation of individual Skolem sequences 
of large order, and the strict enumeration of Skolem sequences. 

Skolem proved in 1957 that Skolem sequences only exist when n = 0, 1 mod 4 [22 . There are many 
methods of generalizing Skolem sequences. These include adding a hook, or 0, into one or multiple positions 
in the sequence, allowing A > 2 copies of an integer, or stipulating that each integer i £ {1, . . . ,n} must be 
present times in the sequence. In |T5], Navena Francetic and Eric Mendelsohn gave a thorough survey of 
the known existence results regarding many of these generalizations. 

Eldin et al. |14j gave a hill climbing algorithm to generate Skolem sequences of arbitrary order, easily 
constructing Skolem sequences of order 84 in examples. 

The enumeration of Skolem sequences has proven quite challenging. In 1986, Jaromi'r Abrham determined 
a lower bound for the number of Skolem sequences using a construction called additive permutations. He 
showed that if n = 0, 1 mod 4, then |iS„| > 2LtJ 1 . Other than such bounds, exact enumeration of the 
number of Skolem sequences has seen most success with inclusion-exclusion algorithms [16 , [24J. 

The exact number of Skolem sequences of order n has been computed for n < 23 only (OEIS A004075), 
a limited number that is mainly due to the very large search space. Before the inclusion-exclusion algo- 
rithms were introduced |16j |24j , enumerative results were restricted to the exhaustive generation of all (2n) ! 
permutations of {1,2,..., 2n}, and counting the valid ones. 

Strategies for exhaustive generation of Skolem sequences have also been quite limited. While [16] and 
[24] used an inclusion-exclusion algorithm to enumerate Skolem sequences of order up to 23, these methods 
did not create the sequences, only counted them. It has been classically believed that the only way of 
generating all possible Skolem sequences of order n is to check all (2n)! permutations. Since, as stated, 
Skolem sequences can be applied to many other areas, the actual construction of all Skolem sequences is very 
useful. The strategy we propose here can theoretically generate all Skolem sequences up to computational 
limitations. In our computing environment, we were able to actually generate all of the sequences up to size 
7, and some of size 8. The main characteristic of this strategy is a significantly reduced search space. 
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4. Skolem sequences as arc diagrams 



As suggested by the original problem posed by Skolem, the partitioning of the set {1,2,..., 2n} provides 
a natural way to depict these sequences using arc diagrams. Vertices are labeled from 1 to 2n of the 2n 
positions in the Skolem sequence and if entries s,; — Sj in w £ <S„, then an arc {i,j) is drawn between vertices 
i and j in the set partition arc diagram. In an arc diagram we define the length of a given arc {i,j), where 
|j — i| = k, as k. Notice that the length of an arc in a Skolem sequence's arc diagram is Si = Sj ~ k. 

Example 3. The Skolem sequence w — (3, 4, 2, 3, 2, 4, 1, 1) G S4 is depicted as the following arc diagram: 




1 2 3 4 5 6 7 8 
w = 3423241 1 

Once a Skolem sequence is represented this way, we can easily notice that it is not only a special case of a 
set partition, but also a special case of a perfect matching. In [3] a generating tree strategy was employed to 
enumerate fc-nonnesting set partitions using a carefully chosen label; here a similar method will be employed 
using open Skolem sequences. However, instead of a label that keeps track of a nesting index, we will track 
arc lengths. 

Definition 2. An open Skolem sequence of order n, p ~ (si, . . . , s„) is a decorated integer sequence made 
up of integers < n, some possibly decorated with a *, such that the following conditions hold: 

(1) If Si and Sj are not decorated with a *, then there are exactly two elements Si and Sj such that 
Si = Sj = k, and \j — i\ = k. In the arc diagram, an arc of length k is drawn between vertices in 
positions i and j. Also, such a k is unique. 

(2) If Si and Sj are both decorated with a *, then Si ^ Sj for j ^i. If Si = *k, then in the arc diagram, 
vertex i — (n + 1) — fc is an open arc that is not closed. Notice that i is the minimum length that the 
open arc must be when completed. 

We denote the set of all open Skolem sequences as OS and those of order n with OSn . Notice that for 
p G OSn, \p\ = n, while for w G iS„, \w\ = 2n, so the corresponding arc diagrams of a open Skolem sequence 
and Skolem sequence, each of order n are n and 2n respectively. 

To use the language of |5, open Skolem sequences can be represented with open arc diagrams, where the 
correspondence between elements in the open Skolem sequence and vertices in the diagram is as follows: 



Vertex type 


Arc Diagram 


Si = 




{ 




semi-opener 


i 


*k 




r . 




opener 


i j 


k 








closer 


j i 


k 



Example 4. Consider the open Skolem sequence p = (*7, 4, 1, 1, *3, 4, *1) £ OS-[ and its arc diagram 
representation: 




From an open Skolem sequence of size n we can construct its descendants through the addition of these 
two distinct vertex types to its arc diagram: 
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(1) An opener (see vertices in position 1,2,3,5 and 7 in Example U]) may always be added in position 

71+1. 

(2) A closer, (see vertices in positions 4 and 6 in Example U]) may be added in position n + 1 provided 
an available semi-opener exists. 

Example 5. Consider two descendants of the open Skolem sequence w G OS-j seen in Example^ that will 
each have order 8: 




12345678 12345678 
7411 *4 4 *2 7 *8 41134 *2 3 



Notice that any descendant of this open Skolem sequence will always have S2 ~ s^ — A and 3^ — 34^ — 1. 

5. A GENERATING TREE FOR OPEN SkOLEM SEQUENCES 

We present a succession rule for the construction of all open Skolem sequences. To each open Skolem 
sequence p G OSn, we associate two labels: 

(1) The p itself; 

(2) A set S made up of the unstarred elements p (without multiplicity). 

Example 6. Take p G OS7 as seen in Example\^ Then p would have labels (*7, 4, 1, 1, *3, 4, *1) and 
5 -{1,4}. 

Skolem sequences may be recognized as the subset of open Skolem sequences of length 2n where S — 
{1, 2, . . . , n}. This may be tested using cardinalities and maximal elements: 

Proposition 1. A Skolem sequence represented with labels of p £ OSn o-nd set S satisfies the following: 

(1) 2|5| = \n\; 

(2) max(S') = IS*]. 

Proof. A Skolem sequence is made up of the elements of the multi-set {1, 1, 2, 2 . . . , n, n}. When k is present 
twice in p G OSn, no * decorates it, and k is in the set S, thus Claim 1 is shown. Since each element 
of a Skolem sequence is present twice, the corresponding label set is — {1,2, .. .n} whose cardinality is 
n — max(S'), thus Claim 2. □ 

We use these labels to generate all open Skolem sequences, and the above proposition to identify Skolem 
sequences. 

We use open Skolem sequences of order n to build open Skolem sequences of order n + 1. There are two 
ways of adding the next element: 

(1) Adding an element with a * decorating it at the end. 

(2) Adding a copy of one of the * elements to alter the * element to a non-decorated state, provided 
that its value is not already present in the S label set: 

In the language of [3 , (1) corresponds to the addition of an opener and (2) to the addition of a closer. 

Example 7. In Examples and ^ we saw an open Skolem sequence of order 7 and two of its children. In 
Example\^ we established that the labels of the given p would be: (*7,4, 1, 1, *3,4, *1) and {1,4}. The labels 
of the first descendant seen in Example^ is (7, 4, 1, 1, *4, 4, *2, 7) and {1, 4, 7}, while the labels of the second 
are (*8, 4, 1, 1, 3, 4, *2, 3) and {1,3,4}. 

We describe the succession rule for constructing open Skolem sequences of order n -I- 1 from an open 
Skolem sequence of order n: 

(1) To add an opener: To each *k £ p £ OSn ~^ *{k + 1), and *1 is appended to p to create a 

p G OSn+l- 

(2) To add a closer: For each *j E p £ OSn, check if j G S. If yes, stop. If no: 

• For all *k ^ *j, *k — > *{k + 1), 
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• *J J, 

• j is appended to p G OSn to create a /i G OSn+i, and 
We depict the start of the generating tree in Figure [TJ 



0, {} 
i 

r 



ri), {} 



(*2,*1),{} 



r r r r ^ .-rs . . r 

(*3, *2, {} (*3, 1, 1), {1} (2, *2, 2), {2} (1, L 1*), {1} 



Figure 1 . The start of the generating tree. 



Example 8. We consider the open Skolem sequences that arise from adding a vertex to p = (5*, 4*, 1, 1, 1*) G 
0^5 in Figure\M 



1 2 3 4 5 
(*5,*4,1,1,*1),{1} 



r r ^ r c 

12 3 4 5 6 
(*6,*5,1,1,*2,*1),{1} 




1 2 3 4 5 6 
(*6,4,1,1,*2,4),{1,4} 



1 2 3 4 5 6 
(5,*5,1,1,*2,5),{1,5} 



Figure 2. w = (5*, 4*, 1, 1, 1*) g OS^ and its descendants. 



6. Discussion 

Programming: The succession rule has been successfully programmed, see or Appendix 1, in Maple 
16 using a 2 Intel Xeon. The generating tree produces open Skolem sequences which are then searched for 
Skolem sequences, using Proposition [TJ The total number of open Skolem sequences is: 

1, 2, 4, 8, 20, 52, 146, 430, 1306, 4176, 13832, 47452, 169044, 619672, 234225 

Note that the sequence made up of |OiSn| is not found in the Online Encycolpedia of Integer Sequences. 
Enumerating the number of open Skolem sequences requires the use of more memory than is given with our 
large computing power, which runs out after over 30GBs are used, and time is less than 10 minutes. Due to the 
nature of our labelling system insufficient memory to compute OSiq is unsurprising, if disappointing. That 
said, while we have not pushed the enumeration of Skolem sequences further, we have dramatically reduced 
the potential search space for future endeavors. For example, historically when exhaustively searching for all 
Skolem sequences of order n = 5, all (2 • 5)! = 3, 628, 800 permutations were considered. Using our method, 
only 4176 open Skolem sequences must be searched to find Skolem sequences of order n = 5. We anticipate 
that there are less than 10 million open Skolem sequences to search in order to find the next exhaustive 
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set of all 504 Skolem sequences of order 8, while classically there would be over 20 x 10"'^^ permutations to 
consider. 

Beyond this purely enumerative benefit of our method, it should be noted that through our generating tree, 
Skolem sequences are also exhaustively constructed, a definite goal in the field of design theory. Furthermore, 
this method lends itself to reasonable modifications that the authors are confident will lead to exhaustive 
generation and construction of other Skolem-type sequences. 

Indeed there are a variety of methods for extending Skolem sequences, including adding in hooks, or O's, 
in the sequence, and allowing more than one pair of entries to have value k. In the latter, this could include 
having entries with the same value k all be k units apart, or different pairs with the same value could be 
independent of each other, see [TU]- In each case, with some careful bookkeeping, we are confident the 
succession rules may be carefully manipulated to include these extensions. 

While the enumeration of Skolem sequences, and various generalizations of Skolem sequences has been 
considered to varying degrees in various algorithms ([10], [24], [16]), the strategy presented here is a departure 
from them in two different ways. Firstly, Skolem sequences of size n + 1 are built from open Skolem sequences 
of size n, and all may exhaustively be generated in this manner. Secondly, and perhaps more importantly, 
this method in which the sequence is represented as an arc diagram and then a label is used in order to 
produce a succession rule for a generating tree speaks to a larger picture in which arc diagrams provide a 
unifying theory. By encoding Skolem sequences in arc diagrams, we may use this framework for construction 
and enumeration. This framework highlights the potential for many other combinatorial classes to be studied 
in this manner. 

Remark: In [5], matchings are encoded in oscillating tableaux with integer fillings in order to prove 
equidistribution between crossing and nesting statistics. A Skolem sequence, when viewed as a subset of 
matchings, corresponds to those oscillating tableaux in which each integer filling appears k times for every 
k e {1,2, ...n}. 

Example 9. We have seen the Skolem sequence w = (3, 4, 2, 3, 2, 4, 1, 1) e ^4 drawn as an arc diagram: 



1 2 3 4 5 6 
3 4 2 3 2 4 



Its filled oscillating tableaux is as follows: 



2 



1 


2 


3 





7 8 
1 1 



We see that the 7 appears once, the 3 appears twice, the 1 appears three times, and the 2 four times. The 
presence ofn different integers, appearing each of the {1, . . . ,n} times is unsurprising for a Skolem sequence. 
This is because an arc of length k enters the tableaux at the presence of its closer, (when read right to left) 
and does not leave until its opener has been reached. 

Future work includes the construction of generalized Skolem-type sequences using this method and the 
translation of the succession rule to a functional equation for faster enumeration. 
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Appendix A. Maple code 

### Preamble ########### 

'type/openSkolemlabel' := [ list (integer) , set(posint) ]; 
incrementstar := procCn: : integer) 
if n >= then 

n; 
else 

n - 1; 

fi; 
end; 

addopener := proc (L : : openSkolemlabel) 

[ [ op (map (incrementstar , L[l])), -1], L [2] ]; 
end; 

addcloser := proc (L :: openSkolemlabel) 

seq( addcl(L, i) , i=l . .nops (L [1] ) ); 
end; 

addcl := proc (L :: openSkolemlabel , i::posint) 
local K; 

if L[l][i] >= or member (-L [1] [i] , L [2] ) then 
return NULL; 
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fi; 

K := L[l] ; 
K[i] := -K[i] ; 

[ [ op (map (incrementstar , K)), -L[l][i] ], 
L[2] union { -L[l] [i] } ] ; 

end; 

Skolem := proc (L : : openSkolemlabel) 
local n; 

n := nops(L[l] )/2; 
type(n, integer) 
cind 

L[2] = {$1. .n}; 

end; 

### end preamble ############### 
L :=[[], O ]; 

to 20 do 

L := opCmapCaddopener , [L])), op(map(addcloser , [L])); 
K := select (Skolem, [L] ) ; 
n : = nops (K) ; 
if n > then 

print ( n, op(K) ) ; 

fi; 

od: 
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